1 

Joint Estimation of Information and Distributed 
Link- Scheduling in Wireless Networks: Mean-Field 
Approximation and Graphical Models 

Sung-eok Jeon, and Chuanyi Ji 
sujeon@microsoft.com, jic@ece.gatech.edu 



Abstract — In a large multi-hop wireless network, nodes are 
preferable to make distributed link-scheduling decisions with 
information exchange only among a small number of neighbors. 
However, for a slowly-decaying channel and densely-populated 
interfering nodes, a small size neighborhood often results in 
nontrivial link outages and is thus insufficient for making optimal 
scheduling decisions. A question arises how to incorporate the in- 
formation outside a neighborhood in distributed link-scheduling. 
In this work, we develop joint approximation of information 
and distributed link scheduling. We first apply machine learning 
approaches to model distributed link-scheduling with complete 
information. We then characterize the information outside a 
neighborhood, i.e., the residual interference, as an aggregated 
random loss variable. The loss variable is characterized by either 
a Mean-Field approximation or a normal distribution based on 
the Lyapunov central limit theorem. The approximated informa- 
tion outside a neighborhood is incorporated in a factor graph. 
This results in joint approximation of information and distributed 
link-scheduling in an iterative fashion. Link-scheduling decisions 
are first made at each individual node based on the approximated 
loss variables. Loss variables are then updated and used for 
next link-scheduling decisions. The algorithm repeats between 
these two phases until convergence. Interactive iterations among 
these variables are implemented by a message-passing algorithm 
over a factor graph. Simulation results show that using learned 
information outside a neighborhood jointly with distributed link- 
scheduling reduces the outage probability close to zero even for 
a small neighborhood. 

Index Terms — Distributed Link-Scheduling, Wireless Net- 
works, Mean-field approximation, Lyapunov central limit the- 
orem, Message Passing, Factor Graph. 

I. Introduction 

Efficient link-scheduling and channel assignment are fun- 
damental to optimal resource utilization in wireless networks 
|5||14||22|. Distributed scheduling is a frequent scenario, 
where individual nodes make scheduling decisions using only 
information from their neighbors [2|[21 1[22|[33|. 

Such distributed scheduling is optimal if an objective func- 
tion to be optimized is local by nature, i.e., decomposable to 
factors with nodal variables and their neighbors 12111221 1331 . 
Distributed scheduling with local information-exchange be- 
comes non-optimal if information outside neighborhoods is 
pertinent. For example, interference causes long-range spatial 
dependence in multi-hop wireless networks [16] [28]; thus op- 
timal scheduling decisions made by nodes require information 
not only within neighborhoods but also outside. 

In our prior work, we studied such a scenario where there 
was a long range spatial dependence among nodes but the 



dependence was weak. We developed distributed scheduling 
based on a machine learning framework, where information 
exchange occurred only among neighbors. Information, i.e., on 
residual interference, was entirely ignored outside a neighbor- 
hood range. We derived a sufficient condition to show that the 
aggregated spatial dependence outside a neighborhood needed 
to be weak, i.e., decay sufficiently fast, for such a distributed 
algorithm to be near-optimal. 

This work studies the scenarios that the spatial depen- 
dence is of long-range and strong. The information, if ig- 
nored outside neighborhoods, would result in a significant 
link outage probability. Therefore, distributed scheduling with 
information exchange only from neighbors is no longer near- 
optimal lfTBI . This corresponds to the scenario of either a 
slowly-decaying channel or densely populated interferers. 
There, the aggregated residual interference cannot be ignored 
|1||3]|7||28||16|. Hence, a node needs to possess informa- 
tion outside its neighborhood, e.g., on residual interference, 
when making cooperative scheduling-decisions. Otherwise, 
decisions on channel-access made by a node may cause 
significant interference even at faraway nodes in a network. 
In our prior work |16|, we show that residual interference 
outside a neighborhood cannot be ignored especially for slowly 
decaying channels and densely populated nodes. However, 
how to incorporate the information outside a neighborhood 
is not studied there. Overall, systematic study is also lacking 
on what size neighborhood is sufficient and what information 
to include for optimal distributed link-scheduling. 

One scenario is for nodes to simply include complete in- 
formation outside their neighborhood. For example, a protocol 
model based on contention lfl2ll assumes that each node main- 
tains the configuration information from an entire network. 
This, however, does not reflect a realistic communication envi- 
ronment ifTBI . where complete information is often unavailable 
at individual nodes. Even if the complete information were 
available, a node would need to exchange information with all 
others [7 ) [9 1 [27 ], which can be prohibitive in a large network. 
As an alternative, we consider aggregated information outside 
the neighborhoods as an approximation. 

The approximation does not consider information from 
individual nodes outside a neighborhood but their aggregated 
contribution. Such approximation needs to satisfy a certain 
optimality in regard to an underlying exact model. Hence, we 
first provide a network model for distributed link-scheduling. 
The model is developed in our prior work ITBI . and based 



2 



on Boltzmann distribution and probabilistic graphical models 
(see [ 10] [ 17] and reference therein). The model includes long- 
range spatial dependence among all nodes in a network, and 
is exact under given assumptions. 

We then consider two approximations to the long-range 
dependence, i.e., residual interference, outside a neighborhood 
range. Both approximations regard the information outside a 
neighborhood as one aggregated variable. In particular, the first 
approximation views the aggregated residual interference as a 
deterministic mean-field. The mean-field parameters optimize 
a variational bound of the probabilistic graph of our network 
model |20|. The parameter values are obtained through solving 
coupled mean-field equations 11771 ll20l . Mean-field approx- 
imations have been used widely in machine learning and 
image processing [17] in the context of graphical models. In 
networking, Mean-Field approximation has been applied to 
loss networks (TBI . These forms of Mean-Field approximations 
are derived based on an underlying model of networks, which 
motivate our work here. For wireless network, Mean-Field 
has been considered as a limiting form of a large number 
of mobile users |4|. This is related to the idea of aggregating 
dependencies outside a neighborhood in this work. Mean-Field 
theory is also used to derive effective link capacity of wireless 
networks 1 19 1. 

The second approximation treats the aggregated residual 
interference as a random variable, referred to as loss vari- 
able. We show that the loss variable can be approximated 
by a normal distribution through the Lyapunov central limit 
theorem (CLT) H321 . The mean and the variance of the 
normal distribution can be learned using information within a 
neighborhood range. The approximations are combined with 
iterative and distributed link-scheduling decisions. In particu- 
lar, a message-passing algorithm on the factor graph ||231 is 
combined with the approximations of the long-range residual 
interference for the distributed and statistical link-scheduling. 
The mean-field parameters and the normal residual random 
variable form particular local functions in the factor graph. 
The Mean Field approximation or the distribution of the loss 
variables are learned iteratively from the available information 
within the neighborhood range. Distributed decisions are then 
made through alternating between approximation and message 
passing over the factor graph. This process repeats until 
convergence. 

We then conduct simulation on network-wide link outage 
probability to complement the analytical derivation of the 
approximation and the algorithm. Simulation results show that 
distributed link-scheduling with an arbitrarily small neighbor- 
hood and aggregated information outside can satisfy the global 
SINR requirements while maximizing the channel reuse. 

The rest of the paper is composed as follows. Section II 
provides motivating examples on distributed scheduling using 
complete, partial and approximated information. In Section 
Hill we provide a problem formulation. Section [IV] provides 
an accurate probabilistic model of distributed link-scheduling 
decisions at individual nodes. In Section [V] we obtain two 
approximations on residual interference, and the resulting 
simplified probabilistic graphical models. In Section I VII we 
derive distributed message-passing algorithm coupled with the 



approximation. In Section lVTTI we provide simulations to show 
the performance of the algorithms using two approximations. 
Section IVIIII concludes this study and discusses open ques- 
tions. 

II. Distributed Link-Scheduling with Loss 
Variables 

In this section, we provide motivation on the importance 
of information outside a neighborhood. We first describe the 
setting of distributed link-scheduling. We then present an 
example of distributed scheduling with complete information. 
We provide an example how the performance degrades for 
distributed scheduling with partial information. This motivates 
two approximation schemes. We provide an overview on the 
approximation schemes to motivate the subsequent sections. 

A. Distributed Link-Scheduling Decisions 

In a distributed setting, a node determines its channel access 
schedule, i.e., whether to access a channel at a time slot, using 
the information within its neighborhood. The information is on 
the "state that includes the relative positions, the transmission 
statues (who is transmitting to whom), transmission power, and 
scheduling decisions of links. For simplicity, this work does 
not consider power control and assumes fixed transmission 
power for each node. 

The decision made is cooperative and based on two criteria. 
The first is whether the SINR requirements are satisfied at the 
neighboring links if a node chooses to use a channel. The 
second is to maximize the spatial channel-reuse which is the 
number of concurrently active links in the network. 

Decisions are made in a distributed fashion. A node sends 
its decision to neighbors. The neighbors use the information to 
make their own link-scheduling decisions, and then send the 
decisions to their neighbors. Distributed decisions are made 
iteratively and collectively by all nodes in a network until 
convergence. 

An assumption is that a node knows complete information 
from its neighbors. Hence, the neighborhood size relative to 
the entire network determines whether the information used 
by a node is complete or partial. 

B. Complete Information 

First, we consider an ideal setting, where complete informa- 
tion is available for individual nodes to make link-scheduling 
decisions |8|[30|. The complete information is on the state 
from all nodes in a network. That is, the neighborhood for any 
node is an entire network. A node then determines whether to 
access a channel at a timeslot using the available information. 
Such a scenario is used as a baseline for comparing approxi- 
mations. 

Consider an example of a linear network in Figure Q] 
Assume that the channel-contention constraint requires any 
active link between two adjacent nodes to be separated by 
one neighboring link. The interference constraint requires any 
active link to be separated by at least two links. Assume 
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Fig. 1. A Linear Topology with 10 Nodes and 9 Directional Links. 

also that the network achieves the spatial channel-reuse max- 
imization where the total number of concurrent active links is 
maximized. 

Figure |2] shows all possible configurations of active links. 
Each row corresponds to a snapshot of active links at a time 
epoch, where the active links satisfy the aforementioned con- 
tention and interference constraints. For example, the first row 
denotes a configuration where both link (3,4) and (8,9) are 
active. Overall, there are 2 configurations with two active links. 
The remaining 8 configurations have three active links that 
result in the maximum channel reuse. Note that for distributed 
decisions where each node has complete information on the 
activities of the others, the configurations with three active 
links are ones that maximize the channel reuse while satisfying 
the constraints. 

C. Partial Information 

Partial (state) information is available to a node when its 
neighborhood consists of only nodes from parts of a network. 

(1,2) (2,3) (3,4) (4,5) (5,6) (6,7) (7,8) (8,9) (9,10) 



a chosen neighborhood size. For example, if two links (2, 3) 
and (7, 8) are active at the beginning, a configuration with 
3 active links can be obtained if the distributed algorithms 
are randomized, i.e., nodes make probabilistic decisions for 
transmission. Our prior work also derives sufficient conditions 
for such distributed scheduling to be near-optimal. However, 
our prior work does not study sufficiently when partial infor- 
mation in distributed scheduling is insufficient for an obtaining 
optimal configuration. 

D. Importance of Information Outside Neighborhood 

The above example shows that for a finite and often small 
number of neighbors, the SINR requirement is not guaranteed 
1 1 6 1 1 28 1 . In other words, information within a neighborhood 
can be insufficient for making correct distributed decisions 
when the information outside the neighborhood is neglected. 

To further illustrate the importance of information outside 
a neighborhood, we consider another example of a multi- 
hop wireless network. The network has nodes located in 
a square area of 10 x 10 square meters, L = 200 links, 
a — 4 as the channel attenuation factor, and SINR threshold 
SINR t /, = 10. Here a = 3 corresponds to slow power 
decay of the channel, resulting in strong interference and thus 
dependence among nodes. Such a scenario is not studied in 
our prior work lfl6ll . where the strong dependence violates the 
sufficient condition we derive for distributed scheduling with 
a moderate neighborhood to be near-optimal. 

Suppose we ignore the residual interference outside a neigh- 
borhood range, which is common practice for large wireless 
networks. We then calculate the link outage probability that the 
SINR requirement is violated at a link. The outage probability 
is shown to be significant and not negligible as in Figure [5] 
especially for a moderate neighborhood size. 



Fig. 2. Possible configurations. 



Active link 



For example, assume the neighborhood for each node is 
the adjacent link. This means that the neighborhood is only 
determined by the contention constraint where a node would 
not access the channel if its neighbor is transmitting. Such a 
neighborhood size is smaller than desirable since node does 
not know whether SINR is satisfied at its neighbors' neighbors. 
Hence, the SINR constraint may be violated. 

Specifically, under the above assumption on the neighbor- 
hood size, it is legitimate for node 3 to transmit to node 4 when 
links (1,2), (6,7) and (9,10) are active (see the eighth row 
in Figure [2j. Those active links are outside the neighborhood 
system and thus invisible to node 2. Active link (3,4) thus 
generates interference, violating the SINR constraints at the 
existing active links. 

By increasing the neighborhood size, we show, in our prior 
work [16], that an optimal configuration with three active links 
can be obtained through fully distributed node decisions with 




Neighborhood range y, 
Fig. 3. Link Outage probability. 

Note that the link outage due to the ignored residual 
interference is different from the transient outage due to the 
fading channel. The link outage here is persistent. Thus, the 
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network connectivity may not be guaranteed. Some nodes can 
even be permanently disconnected if outages occur persistently 
at adjacent links. 

This example clearly shows that the residual interference 
cannot be arbitrarily ignored in a large wireless network. 



E. Approximation of Information Outside Neighborhood 

Now consider that each node accounts for information 
outside its neighborhood, i.e., the residual interference, as a 
loss variable. The loss variable treats the impact from the rest 
of a network except the neighbors in aggregation. 

Two approaches are used for the approximation. The first is 
to characterize the expected value of a loss variable through 
the deterministic Mean-Field approximation 11241 . That is, a 
set of Mean-Field equations are obtained through optimizing a 
variational bound from the network model [20|. The expected 
values of the loss variables can be obtained iteratively through 
solving deterministic Mean-Field equations. The second is to 
regard the loss variable as random and estimate its probability 
distribution. The probability distribution can be learned itera- 
tively using information within the neighborhood. 

The approximations are incorporated into distributed 
scheduling decisions. Each node first estimates either the 
Mean-Field value or the distribution of the loss variable. 
The node then determines channel-access schedule using both 
complete information within its neighborhood and estimated 
information on the loss variable. The algorithm is to be 
described in Section [VT] 



Rij is the inverse of SINR of an active link at a time 
instance, 



R; 



v.. 



(1) 



PiG(i,j) 

where rjj denotes the noise power at the receiver j. For 
simplicity, all nodes assume to have the same (random) noise 
power rjj = To satisfy the SINR constraint, Ry < R t h 
holds, where R t /j = l/SINR t /j is the inverse of a required 
SINR threshold SINR t /j. Let Ry be the approximation of Ry 
that includes complete information within the neighborhood 
of node i plus the residual interference outside, 

rjj + ReSjj 



E|x m -X3'|< 7 f P m G(m,j)a r , 
R, > — 



PiG(i,j) 



(2) 



where ReSy denotes the approximated value of residual in- 
terference outside the neighborhood range experienced by the 
receiver j of active link (i, j). The distance between the closest 
interferer and the receiver of the active link is referred to as 
the contention range r c of an active link. 

B. Objectives 

The objective of link scheduling is to maximize the spatial 
channel-reuse at a time instance, while satisfying the SINR 
requirement of active links, i.e., 



maximize 



subject to 



(Ro - Ry)crj 



E 

Rth >R 



(3) 



y , for V cry- 



Ill. Problem Formulation 

We now obtain a problem formulation for distributed 
scheduling with loss variables. 

A. Notations 

Consider a wireless multi-hop network with n nodes. Xi and 
Pi denote the position and the transmission power of node i, 
respectively. Let X = {Xi, ■ ■ ■ ,X n }, where bold characters 
are used to denote vectors in the rest of the paper. Let L denote 
a set of communication links that access a common wireless 
channel. Let <7y represent the activity of link G L at 

a time instance: o~ij=l if link is active, i.e., using the 
common wireless channel; <7 y =0, otherwise. Let a = {cry} 
where G L. Here node i is assumed to be a transmitter 
whereas node j is a receiver. 

We assume that each node maintains complete information 
within a neighborhood range. The shape of the neighborhood 
for each node can be arbitrary rather than circular. The neigh- 
borhood range can be different for all nodes but bounded by 
a radius 7/. Hence, without loss of generality, 7/ is assumed 
to be the neighborhood range for each node. 

The channel gain for link between node i and j is 
denoted as G(i,j). For simplicity, G(i,j) — \Xi — Xj\~ a , 
where a is the channel power attenuation factor and 2 < a < 
6. G(i, j) is random when positions Xi's are assumed random. 



where Ry is the inverse of the interference at receiver j from 
(H). Rth is the inverse of an interference threshold. Rq is a 
large positive constant where Rq > Ry for V G L. Rq 
is used for the objective function to be maximized whenever 
a new link becomes active. The first term £\ . ay is the total 
number of active links which should be maximized for the 
channel reuse in the network. The second term Ry Uy 
represents the effect of total interference from active links. 

The objective function in ([3]) satisfies the following three 
properties. The objective function is (i) constrained to satisfy 
the SINR requirement of any active links; (ii) monotonically 
increasing with respect to the number of active links; and (Hi) 
monotonically increasing as the inverse SINR of an active link 
decreases. Therefore, the two terms in the objective function 
make trade-offs between the channel reuse and the SINR 
requirements. 

To measure the optimality of the objective function, we use 
the link outage probability P(SINRy < SLNR^). As soon to 
be shown, a limited neighborhood size is a major cause of 
significant link outage probability considered in this work. 

C. Probabilistic Decisions 

Optimization is achieved through distributed decisions at 
individual nodes. Nodes determine whether and when cry 's 
become active. Decisions made by a node need to account 
for loss variables, i.e., the aggregated information outside its 
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neighborhood, and decisions made by neighbors. As deter- 
ministic distributed decisions can be easily trapped in local 
minima, we consider probabilistic, i.e., randomized decisions 
(see ifTUl lfl3l lfl6l and references therein). At each time 
epoch, nodes make the link-scheduling decisions based on a 
probabilistic model of the neighboring nodes' link-scheduling 
decisions at the previous time epoch and the loss variable. This 
requires a probabilistic decision model. 

IV. Probabilistic Model 

The probabilistic model of distributed link-scheduling deci- 
sions is based on a machine learning framework in form of a 
Boltzmann distribution and a dependency graph. The model 
is drawn from our prior work [16], and included here for 
completeness. 

A. Boltzmann Distribution 

Our model considers ov, 's as a random field of link activ- 
ities. Such a random field resembles a particle system with 
interacting particles (spins) [34|. States of particles can be 
viewed as binary random variables. A configuration lu of 
interacting binary-state particles can be described by a system 
potential energy H(u>), The probability distribution of H(ui) 
obeys the Boltzmann distribution P(oj) = Z^ 1 - exp(— H(oj)) 
lfT0lll34l . where Zq = ^ w exp (— H(ui)) is a normalizing 
constant. 

We now regard the objective (01 as the system potential 
ener gy and appeal to Boltzmann law in statistical physics 
[ 10 1 [16 1 1 25 1 1 33 1 [34 1 . In particular, the system potential energy 
of binary link-scheduling decisions is defined as, 

H{<t) = ^(- J Ro+R lJ )a 1J +,3^£/(R lJ -Rthha (4) 



A desirable configuration corresponds to a small system 
potential energy, and a large Boltzmann probability iflOll . 
Hence, an optimal decision cr* is the one that maximizes the 
Boltzmann distribution, 



cr* = argmaxP(er). 



(6) 



B. Dependency Graph 

The probability model exhibits a graphical representation 
1 1 6 1 1 17l ll23l . The graphical representation shows the statistical 
spatial dependence among link activities. For ease of illus- 
tration, consider a directional linear topology of eight nodes 
and seven directional links, where the directions correspond 
to a data flow. The corresponding dependency graph of link- 
scheduling decisions cr is shown in Figure |4] 

Nodes in the graph represent binary random variable {cfij} 
for cry 6 {0, 1} and 1 < i < 7 and j = i + 1. A link 
between two nodes represents the spatial dependence. The 
spatial dependence results from interference in the wireless 
channel and corresponds to coupled terms in the system 
potential energy H(cr). The dependency graph can be further 
represented in form of a factor graph shown in Figure [5] A 
factor graph consists of variable nodes and functional nodes 
ll23l . Here, cr^ denotes a variable node for link S L. A 
functional node tf> i . = Rij denotes a potential function of ov, . 



+ 



E E 



p m i: 



Pii~° 



+ E 777^ + p E u ^ - 



Pd 

13 - 13 

where j3 > is a constant. f3U(Rij — Rth) is a penalty term 
for the SINR constraint. Nij denotes the set of links within 
the neighborhood range from receiver j of link and hij 
denotes total residual interference outside the neighborhood 
range experienced by active link i.e., hij = J2 mn ^N - 

p r a <7 

m mj Tan- 

From Boltzmann law [10|[16|, the probabilistic model for 
a set of link-scheduling decisions at a time instance is, 



P(ct) = Z- 1 - S k P (-H(ct)), 



(5) 



where cr is a random field (i.e., configuration) of all link 
activities cr^'s. Z a = J2cr ex P(^H(cr)) is a normalizing 
constant. Note that both the number of active links and 
interference play important roles: If multiple configurations 
exist with the same number of active links, the system potential 
energy distinguishes different configurations by the sum of the 
inverse SINR of active links. 



Fig. 4. Dependence Graph of random field cr 



Fig. 5. Example Factor Graph Representation. 

Link-scheduling decisions are all dependent due to interfer- 
ence. Thus the dependency graph is fully connected. Such a 
fully connected graph shows that the neighborhood of a node 
is in fact the entire network. 

C. Optimal Configuration and Distributed Decisions 

A fully connected dependency graph impacts how dis- 
tributed decisions should be made using the probabilistic 
model. Node decisions are all dependent. Hence, optimal 
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decisions require nodes to exchange information with all the 
other nodes in the network. Graphically, this means that an 
entire network is a neighborhood. That is, a node needs to 
obtain information on all other nodes. This, in reality, prohibits 
implementation of optimal distributed link-scheduling. 

V. Approximation of Probabilistic Model 

We now derive approximations of the probabilistic model 
for distributed scheduling. Our approximation considers a 
moderate neighborhood for each node where the node has 
complete information on its neighbors. Outside the neighbor- 
hood, the residual interference, instead of being ignored, is 
characterized as a random loss variable. Either the mean values 
or the distributions of the loss variables are then included 
in the probabilistic model. The resulting model P l (<r) is an 
approximation of the exact model P(cr). 

A. Approximation by Mean Field 

We first characterize the residual interference by Mean 
Field (MF) approximation. The MF approximation optimizes 
a variational bound between the approximated and the exact 
models ll20l . resulting in a set of deterministic mean-field 
equations. The mean-values of the loss variables can then be 
obtained through solving the Mean-Field equations iteratively. 
To derive the Mean-Field equations, we consider the system 
potential energy in (|4j. For simplicity, we assume that the 
SINR constraint is satisfied so that the system potential energy 
is reduced to 



ij ij mn^Nij 



1-3 



where a 



-Ho 



hij — ... v • 



^pa, »«j,m n - pd - a from ©, and 

ij 1 ij 

mngNij »i3,mn<rij<7mn corresponds to the aggregated 
residual interference outside the neighborhood for link 
Note that this system potential energy represents a second- 
order Ising model in statistical physics, where all terms are 
dependent. 

Let h*j be a deterministic approximation of hij. An ap- 
proximated system potential energy absorbs the impact of the 
aggregated residual interference into the coefficient a^'s, 

#V) = (8) 

ij ij mn£_Nij 

where by = a^- + and h*j is yet to be found. 
We denote the true free energy 1201 with F, 

l3F = -ln(Z P ), (9) 

where /3 > is a constant, and Zp = J2cr exp(— 7J(<x)). Let 
variational free energy [20 1 be F, 

pF = -ln(Z Q ), (10) 



where j3 is a large positive constant, and Zq = 

£ exp(-ffV)). 

To find the simplest, i.e., the first-order Mean Field (MF) 
approximation, we make the following assumptions: (a) for 
link ^ (m,ri), the impacts only affect the first order 

terms; (6) for (m,ri) £ N^, Uij is independent of a mn ; 
and (c) for the approximated system potential energy, the 
corresponding distribution is C](er) = Zq 1 exp(— H l (cr)). The 
mutual information measure between the approximated and 
the true distribution results in an objective function, i.e., a 
variational free energy |20l . 

PF = E Q (ln(^))-HZp) (ID 
= E Q (\n(H(a))-\n(H l (a)))-\n(Z Q ) 



= -HZq) + 



1 f-l'ij f-Ln 



where fiij=E(aij)=P(aij — 1). Inserting Equation (fTTTb in 
(fTTb and letting ^^-=0, the optimal value of hij satisfies 



mn^Nij 



(12) 



h*^ is a weighted sum of expected values of link activities 
outside the neighborhood range. Hence, H*a can be regarded 
as the mean field outside the neighborhood of link As 
a special case when aij, mn = for van Nij, h*^ = 0. 

For £ L, there are n such equations, where n is 

the total number of nodes in a network. Those equations are 
coupled. The coupling can be shown explicitly by rewriting h*j 
(7) as h'lj = log —•■ Equation ( fl2] i then takes the following 



form, 



lot 



1 - fHj rji 



lZ a 

mn^Nij 



ij : mnf^mn j 



(13) 



where ?/i and r\2 are functions of h* kl 's for (fc, I) ^ 
(and in turn phi's)- Hence, the set of n equations are coupled 
through unknown values of h*j's or /i^ 's. AppendixUprovides 
detailed derivations. 

The coupled Mean-Field equations can be solved iteratively 
and numerically starting from an initial condition |20|. The 
Mean-Field entity h*j can also be estimated when measure- 
ments are available. For example, the interference within 
the neighborhood range can be known with the information 
exchange. Thus, if we assume that the aggregated interference 
can be measured possibly with a measurement error, h*j can 
be estimated when interference measurements are available. 

An advantage of the Mean-Field approximation is that 
a resulting probability model is a simple Markov random 
field, where the deterministic Mean-Field values h*ij& absorb 
the long-range spatial dependence in the residual interfer- 
ence. However, the Mean-Field approximation exhibits two 
disadvantages. First, it can be computationally intensive to 
solve a set of Mean-Field equations. Second, if the Mean- 
Field parameters are estimated by measurements instead, such 
measurements on the aggregative residual interference are 
either unavailable or computationally intensive to obtain. For 
example, such measurements require nodes to transmit test 
signals during link-scheduling decisions. In addition, nodes 
need to be capable of measuring the aggregated interference. 
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For example, packet reception ratio (PRR) versus SINR mod- 
els need to be seeded by 0{n) trials in an n-node network 
where each node transmits while receivers measure the channel 
conditions l26ll . 



and r u , for fc ^ m. 

The proof can be found in Appendix HI1 This lemma shows 
that random positions of interferers outside a neighborhood 
follow similar characteristics to within the neighborhood. 



B. Approximation by Lyapunov Central Limit Theorem 

An alternative approach is to treat the aggregated residual 
interference as random and approximate its distribution. Such 
an approximation can be computationally less expensive to 
obtain than the Mean-Field approximation. 

Let riij be a random variable that approximates the residual 
interference of an active link Our goal is to identify 

a probability distribution for n.y. As residual interference 
aggregates interference signals from many interferers, one 
option is to approximate ny as a normal random variables. 

We first obtain an expression for to understand whether a 
normal approximation is feasible. Based on our assumptions, 
interferers have the same transmission power Pq. Locations 
of interferers outside a neighborhood range can be arbi- 
trary/random but follow certain patterns of node-positions 
within the neighborhood. In other words, the network topol- 
ogy exhibits a certain regularity. Specifically, the following 
assumptions are posed on the location of interferers. 

• For an active link all interferers are located uni- 

formly within radial bands bounded by concentric circles 
from receiver j. The radius of the circles exhibits random- 
ness, i.e., r k = Tk~i + x k for k > 1, where ro=7/, and 
-x is a uniform random variable between r l and r". 



Here, r and r" are the minimum and maximum distance 
between any two active links, within the neighborhood 
of link The circle with radius r k is called the fc-th 
frontier. 

On the fc-th frontier for k > 1, any two neighboring 
active links are separated by random distance Xq. Random 
variables x k and xq are independent and identically 
distributed. 



1h . 



Fig. 6. Configuration of Active Links outside Neighborhood range. 

Lemma 1: With the above two assumptions, the residual 
interference hij for active link can be modeled as 



^(P r- Q ) • (2nr k /x ) 



(14) 



fc=i 



where r k = r k -i + x k with ro=7/ as a constant. x k = x is a 
random variable for k > 1. x and xq are both independent and 
identically distributed with a uniform distribution between r l 



Next, we derive the mean and the variance of ny's. For ana- 
lytical feasibility of obtaining the mean value, we approximate 
the joint probability density P(rx, ■ ■ - , Too) with a product 
form based on the Mean-field approximation ll24l . 

Lemma 2: Consider {ri, ■ ■ ■ , r^} as a random field. Based 
on the first-order Mean-field approximation, the joint proba- 
bility density function P(ri, ■ ■ ■ , Too) can be approximated as 
a product of marginal distributions, 



n p (nl<n-i >), 



(15) 



where < rj_i > = 7/ + (£ — l)r c denotes the expected value of 
Ti-\- r c is an expected value of r c , i.e., r c = E(r c ) = r , 
and < r >=7/- 

The proof can be found in Appendix!]]]] Using such a mean- 
field distribution, we obtain the expected value 



ri, ■••,r DO J x fc=1 



(Por k 



x 



(16) 



P{r k \ < r fc _i >)P(x )d ri ■■■d rao d x 



where P{r\\ < r >)=P(ri\^f). Closed-form expressions 
can be derived in the following lemma for the individual terms 
in the sum. 

Lemma 3: Let u k be the term in the summation in ( [ToT l. 
where u k = ^P^r\- a = 27rP (r fe _ 1 + x k f- a . 

Denote the mean and the variance of u k as E{u k ) and 
V(u k ) respectively. For a > 2, 



E(u k ) 
V(u k ) 



2itP a (V\- a - V\- a ) 

(2-a)(r u -r l ) 
E(ul)-E(u k f, 



where T>i = 7/ + (fc — l)r c + r u , T>2 = 7/ + (fc — l)r c + r l 



and E(ul) =4tt 2 P 2 {Vr' M - V^ a )/((3 - 2a)(r u - r 1 )). 
Let the fc-th term of ray be v k = Then, E(v k ) — 

g( "j 1 ° ( ; 7r ' ) . and V(v k ) = 

The proof is given in Appendix II VI Using the results from 
Lemmas 1-3, we obtain the mean and the variance of the 
residual interferences below. 



Theorem 1: For a link residual interference 



is a finite sum of independent random variables with 
different probability distributions. Thus, the probability 
distribution of ny is approximately normal with mean £(ny) 
and variance V(n,ij), where E(riij) — X^feLi E( v k), an d 

The proof is provided in Appendix [V] 

How accurate is the normal approximation? We check the 
modeling accuracy using Lyapunov condition. For ny to 
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converge to a normal distribution, the Lyapunov condition on 



J2T=i v k should satisfy 

(£r =1 £(^-£K)i 3 )) 1/3 



1/2 



(17) 



Such condition is valid for an infinite sum of independent 
random variables. As the interference from far apart interferers 
is diminishing to zero, riy is in reality a sum of a finite number 
of independent random variables. As a result, the Lyapunov 
condition converges to a small constant instead of zero. A 
numerical analysis on the normal approximation is conducted 
on a network with infinite nodes, the increase of neighbor- 
hood size to fc-th frontier neighbors, and channel attenuation 
parameter a = 4. Figure [7] shows that the Lyapunov condition 
converges to a small constant around 0.97. 

Hence, the advantage of the normal approximation is 
the simplicity. Only the mean and the variance need to be 
estimated. However, the approximation can deviate somewhat 
from a normal random variable measured by the Lyapunov 
condition. 
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Fig. 7. Lyapunov Condition of the residual interference rt{j 



C. Approximated Potential Energy and Probabilistic Model 

The Mean-Field and the normal approximations result in 
approximated two network potential energy functions and thus 
the probabilistic models respectively. First, the Mean-Field 
approximation in Section IV- Al characterizes the residual inter- 
ference as a deterministic quantity outside the neighborhood. 
The resulting system potential energy is 



where h* 



(18) 



ij 



is deterministic and obtained from the mean-field 



equation. The corresponding Boltzmann distribution is 

P l (a) = Z- 1 exp(-i/'H). (19) 

Such a Boltzmann distribution is a second-order Markov 
Random Field with deterministic coefficients ay's, h*j's and 



aij,mn's- The resulting probabilistic dependency graph is 
shown in Figure [8] 

Second, the normal approximation characterizes the residual 
interference outside the neighborhood as a Gaussian random 
variable that aggregates contributions from (unseen) interfer- 
ers. The resulting system potential energy is, 



H\a,n) = J2(- R o + 



N h 



P l~ a ' 



u iju rnn 



(20) 



PJ-« aXj 



P l~ a 

^ ^ P,l~, a 

ij mnGN i j LJ 

+/3^[/(4 i -R« s Ki, 



where n is a random vector including all n^'s. The corre- 
sponding Boltzmann distribution is 



P l (cr : n) = Z,7 1 exp(-7j'(cr,n)) 



(21) 



This approximated model is a random-bond Markov Random 
Field lfl6l . where n^'s are random coefficients (bonds) for 

(Jy'S. 

As distributed decisions only need to determine the values 
of CTy's, the marginal probability model for cry's is 



P\rr) = J P l {(T,n)dr 



(22) 



VI. Message-Passing over Factor Graph 

A simplified probabilistic model P l (cr) from either the 
mean-field or the normal approximation can be used to de- 
termine link activities, 



er = argmaxP (er). 



(23) 



Randomized and distributed decisions implement such op- 
timization through message passing over localized factor 
graphs. The general message passing algorithm is provided 
in ||231 and the references therein. The novelty here is to 
couple message passing with approximation of the residual 
interference. The corresponding functional nodes and variable 
nodes are chosen accordingly to couple the approximations 
with distributed scheduling. 




%}) >l'- 
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Fig. 8. Example Factor Graph with Functional Nodes and Variable Nodes. 



9 



A. Message Passing with Normal Residual Interference 

The local potential function in message passing, 
P(i, -\<J N i ) for link (i,j), exhibits a special form from 

approximating residual interference as a normal random 
variable, i.e., 



= N(E(4 j ),V(4 j )), 



(24) 



where ^\ 3 =R\ y N(E{^ 
bution with a mean 



' denotes a normal distri- 



« Si 



r _ . ^\JHi), 



Poi; 



and a variance (from (O) 



fain 



(25) 
(26) 



Message from a function node ip _ to a variable node <Xy is 



P(<7ij = l\ip i3 ,a N i ) = U(R th -ipij), 
where U(R t h - R\j)=h if Rth > R\f 0, otherwise. 



(27) 



Message from a variable node cry to a function node i/r 



is 



^ .) ; ' 



otherwise, 



(28) 

where Q(x)= J Q -j^=e 2 dt. This equation is obtained from 
an equality of P{o l3 = l\a N i.) = J P{o~ij = l\ip l ij} a N r.) 

Consider an example in Figure [8] The second node at the 
top row determines the scheduling of link (2, 3). This node 
maintains only the local factor graph related to link (2, 3). 
The local graph has variables ip l 23 , eri2,<723 and 034. 

The algorithm includes two steps: (a) estimating the param- 
eters for the mean and variance of the residual interference, 
(b) using the estimated parameters in the local functions for 
message passing. For example, in Figure [8] a functional node 
^[2 estimates the mean and variance of residual interference 
ni2, E(n\2) and V(ri\2), from the received message from 
neighboring variable nodes CT23 and 034. The estimation is 
from Theorem 1. 

The functional node ip l 12 sends to a variable node a 12 a 
message P(cr 12 — 1|^4 2 ) a N[ 2 ) = U{R t h — ^12)- The variable 
a 12 sends to all the neighboring functional nodes, ^\ 2 , ^23 
and ^34, a message P(oyi = M^n 1 ) m 



B. Message Passing with Mean-Field Residual Interference 

The local function with the Mean-Field approximation of 
the residual interference, i.e., /i* in ©, is relatively simple, 



P(V>lj = a Q \a N i ,h*j) = 6(ipij - do), 



(29) 



where ^\j=R\j m © with Resy=/i*- in (|S). cio is a constant. 



S(ij)lj - a )=l, if ijjlj = a ; 0, otherwise. 



Message from a function node tj/. to a variable node ov, is 



P((Tij = l\lp i: 



• Kj) = u(Rth - 



(30) 



where U (x) is a unit step function. 
Message from a variable node <xy to a function node ijr 



is 



P{oij = l\a N i ,h*j) = U(Rth — Vii) • 



(31) 



The overall algorithm includes two steps: (a) estimating 
the mean-field parameter h*j's, and (b) obtaining decisions 
(Tjj's through message passing. These two steps alternate until 
convergence. 

VII. Simulation 

We now evaluate network-wide link outage probability to 
assess the benefit of the approximation on link-scheduling. 
As the link-outage probability is not analytically related yet 
with the network model, we conduct simulations for the shared 
channel access on a multi-hop wireless network. 



A. Simulation Setup 

Network nodes are positioned uniformly in a square area 
of 100 square meters, and composed of L = 200 links. We 
consider 2 < a < 6 for the channel attenuation factor. Here, 
2 < a < 3 results in 'strong' dependence among interferers. 
This corresponds to the scenario which is unsolved in our 
prior work. We choose 10 < SINR f /, < 100. We consider the 
maximum spatial channel-reuse at a time instance. 

B. Simulation Results 




1 2 3 4 5 6 7 

Number of Iterations until Convergence 

Fig. 9. Convergence of Distributed Link-Scheduling Decisions 

Using the above experimental setting, we first study through 
simulation whether distributed decisions converge sufficiently 
fast. Figure [9] shows that the distributed link-scheduling deci- 
sions converge through interactions among neighbors in a few 
iterations. 
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We now compare the performance of the following al- 
gorithms: (a) distributed link scheduling with the Mean- 
Field approximation, (b) distributed link scheduling with the 
normal approximation, and (c) the conventional distributed 
link-scheduling algorithm [2|[9][16] that simply ignores the 
residual interference outside neighborhood range. The link 
outage probability is used as a performance measure. Such 
a performance measure assesses the probability that a link 
cannot satisfy the SINR requirement. Here, the link outage 
is mainly caused by the limited neighborhood size. There, 
improved performance signifies the importance of proper 
approximations of information outside a neighborhood. 

Figure [10] plots the outage probability as a function of 
neighborhood size 7/ assuming that all nodes have the same 
size neighborhood. This figure shows that our two approxi- 
mations for the information outside neighborhoods result in a 
small link outage probability given for a small neighborhood 
range 7/. The conventional scheme, however, has a large link- 
outage probability for a small 7/. This shows that for a small 
neighborhood, learning the residual interference with a loss 
variable provides a big gain. The performance with the mean- 
field approximation varies with respect to the measurement 
error of residual interference. When the measurement error is 
small, the Mean-Field approximation results in a small link- 
outage probability across different neighborhood sizes. Over- 
all, the two approximation approaches outperform significantly 
the method that ignores the residual interference, showing 
the importance of approximating the information outside a 
neighborhood. 



It is worth mentioning that for small a, e.g., a — 3, the 
two approximation approaches significantly outperform the 
conventional approach. The scenario of small a is shown in 
our prior work lfl6l to be difficult for distributed scheduling 
with information only from neighbors. Hence, the result here 
shows the importance and the ability of the approximations 
to account for the aggregated residual interference outside a 
neighborhood. 



□ Without Loss Variable 

• With Lyapuanov CLX 

A With MF (Measure Error Var=0.001) 

o With MF (Measure Error Var=0.01) 
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Power Attenuation Factor a 
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Fig. 11. Link Outage probability with 7/=4 and SINR i / l =10. Dotted line: 
The conventional scheme. Other curves: The two approximation schemes. 



□ Without Loss Variable 

• With Lyapuanov CLT 

A With MF (Measure Error Var=0. 001) 

o With MF (Measure Error Var=0. 01) 




Neighborhood range y f 

Fig. 10. Link Outage probability with SINRif^lO. Dotted line: Conventional 
scheme. Other curves: Our approximation approaches 

Figure QT| shows how the link-outage probability varies 
with respect to channel attenuation factor a. The link outage 
probability for the two approximation schemes remains small 
for different a. As a increases, the residual interference 
outside the neighborhood decrease, and the link outage prob- 
ability of the two approximation schemes approaches zero. 



Figure Q~2] and Qj] show the performance as a function of 
SINR. As SINR threshold increases, the link outage probabil- 
ity increases proportionally except for the Mean-Field approxi- 
mation with small measurement error. The performance of the 
Mean-Field approximation for a large measurement error is 
worse as the SINR threshold increases. This is because the 
SINR constraint is more easily violated due to the measure- 
ment error as the SINR threshold increases. 

For the approximation scheme with Lyapunov central limit 
theorem (CLT), once 7/ is moderately large, the proposed 
scheme copes well with a large SINR threshold. When neigh- 
borhood size jf is not sufficiently large, for a moderate SINR 
threshold, conventional schemes in the literature often fail to 
satisfy the SINR requirement, resulting in a high link outage 
probability. On the contrary, across a wide range of SINR 
threshold values, the proposed scheme with Lyapunov CLT 
performs efficiently and improves the link outage probability 
significantly. In particular, if 7/ is moderate (e.g., 7/ = 4 
in Figure [13] ). the link outage probability quickly approaches 
zero. Therefore, with both approximation approaches, resource 
utilization is significantly improved overall. 

VIII. Conclusion 

We have studied in this work distributed scheduling de- 
cisions where each wireless node uses complete information 
within a neighborhood and inferred information outside. We 
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□ Without Loss Variable 

• With Lyapuanov CLT 

A With MF (Measure Error Var=0.001 ) 

o With MF (Measure Error Var=0. 01) 




A— 



20 40 60 

SINR Threshold 



Fig. 12. Link Outage probability with 7/=3. Dotted line: Conventional 
scheme. Other curves: Approximation schemes. 



□ Without Loss Variable 

• With Lyapuanov CLT 

A With MF (Measure Error Var=0.001 

o With MF (Measure Error Var=0.01) 



SINR Threshold 



Fig. 13. Link Outage probability with 7/=4. Dotted line: Conventional 
scheme. Other curves: Approximation schemes. 



have characterized the complex spatial dependence of dis- 
tributed link-scheduling decisions with a factor graph. The 
factor graph enables information exchange among neighbors. 
More importantly, the long-range dependence in residual in- 
terference is not ignored but approximated by one variable 
as the aggregated dependence outside a neighborhood. The 
Mean Field approximation considers the residual interference 
outside a neighborhood as a deterministic Mean-Field. Such 
Mean-Field approximation optimizes the variational free en- 
ergy, resulting in a simple Markov Random Field model for 
distributed decisions. A disadvantage is that the Mean-Field 
estimates can be computationally costly to obtain. The normal 
approximation considers the residual interference outside a 
neighborhood as a random loss variable. Such a loss variable is 



approximately normal by the Lyapunov central limit theorem, 
resulting in a simple Random Bond model for distributed 
scheduling decisions. A disadvantage is that the condition 
for the asymptotic normality can be somewhat deviated for 
a network with a finite size. 

The approximated dependencies are easily incorporated into 
the factor graph. This results in an extended message-passing 
algorithm with the two phases. First, the approximations are 
obtained iteratively from the available information at neigh- 
boring nodes. Second, the approximations are combined with 
messages from neighboring nodes and allow each node to 
make its link-scheduling decisions. The algorithm alternates 
between the two phases, and has been shown to converge in 
a few than 10 iterations. 

We have shown through simulations that the link outage 
probability of the resulting link-scheduling decisions is sig- 
nificantly improved by the approximated information. The 
improvement is especially pronounced for distributed deci- 
sions with a small number of neighbors and a slowly-varying 
channel, where the link-outage probability is reduced from, 
e.g., 0.3 to 0.04. 

An open issue as possible future work is to derive an 
analytical relationship between the network-wide link outage 
probability and the approximated information. 

Appendix 

I. Derivation of Equation (fT"3l 

Hij — P((Jij = 1), where the probability is evaluated using 
the simplified Boltzmann distribution Q(cr) from the mean- 
field approximation. Expanding the expression of /i^ , we have 

exp(-/i*)7?i 



My 



exp(- 



(32) 



where 771 = exp(-ay) £(m,n)^(i,j) exp(-i?(<7 \ o-y)) with 
er \ (Tij being all a mn 's except for (m, n) = rj 2 = 

E( m ,„)/(ij) exp(-if [a\aij = 0)). 

Representing h*j from the above expression, we have 

1 



My Vi 



My V2 

Replacing this expression in ( fT3l l, we have 

1 - My m 



log 



My V2 



E ■ 



(33) 



(34) 



Note that rji and 772 are functions of /i^„'s (MyS) for 
(m,n) ^ (i, j). fiij is also a function of ft^'s. Hence, the 
equation is coupled. In addition, there are multiple equations 
that contain common h*j's. Hence, for all (i,j)'s, the mean- 
field equations are coupled, and can be solved iteratively. 

II. Proof of Lemma 1 

On the fc-th circle for k > 1, the total number of active 
links on the circle is 2zlDl which is random. Each active link 
on the fc-th circle results in as much interference as P r^ a . 
As a result, total interference from the active links on the fc-th 
circle is (P or ^).^; and, thus ny = ££° =1 (P r^ a )-^. 
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III. Proof of Lemma 2 

Mean field approximation shows that the joint distribution 
of a random field (yi, ■ ■ ■ ,y n ) can be approximated by a 
product of marginal distributions. Thus, P{y\ ,•••,?/„) = 11™= 1 
P{Vi\ < Vn % >) with a good approximation, where y^ i 
is the set of neighboring random variables of yi. From the 
assumptions, r fe =r fc _i + x k . Thus, P(r k \ < r Nk >)=P(r k \ < 

n-i >)■ 

IV. Proof of Lemma 3 

The average of u k , denoted with E(u k ), can be derived 
from an integration of J 2itPQr],~ a P{r k \ < r k -i >)dr k . The 
average of u\ can be derived in a similar way. 

V. Proof of Theorem 1 

Consider the configuration in Figure [6] For an active link 
(i, j), the residual interference of this configuration is denoted 
with a random variable mj, which is a function of random 
variables r k for k > 1, where r k is the radius of the fc-th 
frontier. 

Note that r k 's are independent random variables with dif- 
ferent uniform distributions. Lyapunov's central limit theorem 
||32l shows that the summation of a large number of indepen- 
dent random variables (even with different distribution) results 
in a normal distribution, riy is a summation of hundreds and 
thousands of independent random variables (i.e., r k 's), thus 
riij can be approximated by a normal distribution. 

Furthermore, a normal distribution can be completely char- 
acterized by mean and variance. The mean and variance of riy 
is denoted with E(mj) and V(riij), respectively, i.e., E(mj) 
= EfcLi E(v k ), and V{n tJ ) = ~£T=i V(v k ). 

References 

[I] J. G. Andrews. H. Claussen, M. Dohler, S. Rangan, M. C. Reed. 
"Femtocells: Past, Present, and Future," IEEE JSAC, to appear April, 
2012. 

[2] P. Bjorklund, P. Varbrand, and D. Yuan, "Resource Optimization of Spa- 
tial TDMA in Ad Hoc Radio Networks: A Column Generation Approach," 
In Proc. Of IEEE Infocom, April 2003. 

[3] V. Chandrasekhar, J. Andrews, A. Gatherer, "Femtocell Networks: A 
Survey," IEEE Comm. Magazine, pp. 56-67, Sept. 2008. 

[4] A. Chaintreau, J. Le Boudec, N. Ristanovic, "The age of gossip: spatial 
mean field regime" In Proc. Of SIGMETRICS, 2009. 

[5] M. Chiang, "Distributed Network Control Through Sum Product Algo- 
rithm on Graphs," In Proc. of IEEE Globecom, vol. 3, pp. 2395-2399, 
Nov. 2002. 

[6] P. Djukic and S. Valaee, "Link Scheduling for Minimum Delay in Spatial 

Re-Use TDMA" In Proc. Of IEEE Infocom, April 2007. 
[7] T. Elbatt, and A. Ephremides, "Joint Scheduling and Power Control for 

Wireless Ad Hoc Networks," IEEE Trans. Wireless Communications, vol. 

3, pp.74-85, Jan. 2004. 
[8] S. C. Ergen, P. Varaiya, "TDMA scheduling algorithms for wireless sensor 

networks," Journal Wireless Networks, pp. 1435-1440, Vol.16 No. 4, May 

2010. 

[9] S. Gandham, M. Dawande, and R. Prakash, "Link Scheduling in Wireless 
Sensor Networks: Distributed Edge Coloring Revisited," Journal of 
Parallel and Distributed Computing, Vol. 68, Issue 8., pp. 1122-1134, 
Aug. 2008. 

[10] S. Geman, and D. Geman, "Stochastic Relaxation, Gibbs Distributions, 
and the Bayesian Restoration of Images," IEEE Trans. PAMI, vol.6, 
pp.721-741, June 1984. 

[II] A. Ghosh, O.D. Incel, VS. Anil Kumar, and B. Krishnamachari, "Multi- 
Channel Scheduling Algorithms for Fast Aggregated Convergecast in 
Sensor Networks," In Pmc.IEEE Mass, 2009. 



[12] Y. Shi, Y. T. Hou, J. Liu, and S. Kompella, "How to correctly use the 

protocol interference model for multi-hop wireless networks," In Proc. of 

ACM Mobihoc, pp. 239-248, 2009. 
[13] B. Hajek and G. Sasaki, "Link Scheduling in Polynomial Time," IEEE 

Trans. Information Theory, vol.34, pp. 910-917, Sept. 1988. 
[14] B. Han, V. S. Anil Kumar, M. Marathe, S. Parthasarathy, and A. Srini- 

vasan, "Distributed Strategies for Channel Allocation and Scheduling in 

Software-Defined Radio Networks," In Proc. of IEEE Infocom, 2009. 
[15] Ka Hung Hui, Dongning Guo and Randall Berry, "Performance 

Analysis of MAC Protocols in Wireless Line Networks Using Statistical 

Mechanics," In Proc. Of Allerton Conference, 2006. 
[16] S. Jeon, and C. Ji, "Randomized and Distributed Self-Configuration 

of Wireless Networks: Two-Layer Markov Random Fields and Near- 

Optimality," IEEE Trans. Signal Processing, vol. 58, no. 9, pp. 4859-4870, 

Sep., 2010. 

[17] M. Jordan, and Y. Weiss, "Graphical Models: Probabilistic Inference," 
Graphical Models: Probabilistic Inference, 2002. 

[18] F.P. Kelley, "Loss networks," The annals of applied probability, 1999. 

[19] V. Marbukh, "Towards mean field theory of wireless networks," ACM 
SIGMETRICS Performance Evaluation Review, 2003. 

[20] D. J.C. MacKay, "Information Theory, Inference, and Learning Algo- 
rithms," Cambridge University Press, 2003. 

[21] M. Ji and M. Egerstedt, "Distributed Coordination Control of Multiagent 
Systems While Preserving Connectedness," IEEE Trans. Robotics, vol. 23, 
no. 4, pp.693-703, Aug. 2007. 

[22] E. Modiano, D. Shah, and G. Zussman, " Maximizing Throughput in 
Wireless Networks via Gossiping" In Proc. of ACM SIGMETRICS / IFIP 
Performance, Jun. 2006. 

[23] F. Kschischang, B. Frey, H. Loeliger "Factor graphs and the sum- 
product algorithm," IEEE Trans, on Information Theory, vol. 47, no. 2, 
pp. 498-519, Feb. 2001. 

[24] H. Kappen, W. Wiegerinck, "Mean field theory for graphical models,", 
Advanced Mean Field Theory-Theory and Practice, MIT Press, pp. 37-49, 
2001. 

[25] G. Liu, and C. Ji, "Cross-Layer Graphical Models for Resilience of 

All-Optical Networks under Crosstalk Attacks," IEEE JSAC Optical 

Communications and Networking Series, vol. 25, pp. 2-17, 2007. 
[26] S. Liu, G. Xing, H. Zhang, J. Wang, J. Huang, M. Sha, L. Huang, 

"Passive Interference Measurement in Wireless Sensor Networks," In 

Proc. of IEEE ICNP, pp.52-61, Oct., 2010. 
[27] H. Luo, S. Lu, V. Bharghavan, J. Cheng, and G. Zhong, "A Packet 

Scheduling Approach to QoS Support in Multihop Wireless Networks," 

ACM MONET, vol. 9, no. 3, pp. 193-206, June 2004. 
[28] R. Madan, S. Cui, S. Lall, and A. Goldsmith, "Cross-Layer Design for 

Lifetime Maximization in Interference-Limited Wireless Sensor Netwroks " 

In Proc. IEEE Infocom, 2005. 
[29] T. Rappaport, "Wireless Communications: Principles and Practice," 

Prentice Hall. 

[30] I. Rhee, A. Warner, J. Min, L. Xu, "DRAND: Distributed Randomized 

TDMA Scheduling for Wireless Ad Hoc Networks," IEEE Trans. Mobile 

Computing, vol. 8, No. 8, pp. 1384-1396, Oct. 2009. 
[31] P. Santi, R. Maheshwari, G. Resta, S. Das, and D. Blough, " Wireless 

link scheduling under a graded SINR interference model," In Proc. of 

IEEE Infocom, 2009. 
[32] V. A. Statulevicius, "Limit Theorems of Probability Theory, " Springer, 

2000. 

[33] X. Tan, "Self-Organization of Autonomous Swarms via Langevin 
Equation," In Proc. of IEEE CDC, pp. 1435-1440, 2007. 

[34] J. S. Yedidia, and W. T. Freeman, " Understanding Belief Propagation 
and its Generalizations," Exploring Artificial Intelligence in the New 
Millennium, ISBN 1558608117, Chap. 8, pp. 239-236, Jan. 2003. 



